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We study the conductivity of a 3D disordered metal close to antiferromagnetic instability within 
the framework of the spin-fermion model using the diagrammatic technique. We calculate the 
interaction correction 5a(uj,T) to the conductivity, assuming that the latter is dominated by the 
disorder scattering, and the interaction is weak. Although the fermionic scattering rate shows critical 
behaviour on the entire Fermi surface, the interaction correction is dominated by the processes near 
the hot spots, narrow regions of the Fermi-surface corresponding to the strongest spin-fermion 
scattering. Exactly at the critical point 5a oc [mm(w, T)] 3 / 2 . At sufficiently large frequencies uj 
the conductivity is independent of the temperature, and 5a oc (r _1 — icj) -2 , r being the elastic 
scattering time. In a certain intermediate frequency range 5a(uj) oc ia;(r _1 — iu)~ 2 . 
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Transport close to quantum criticality fascinates and 
challenges researches in many fields of condensed mat- 
ter, ranging from physics of high-temperature super- 
conductors and heavy-fermion materials to the conduc- 
tion in graphene or granulated superconductors near the 
superconductor-insulator transition [TJ |2]. At low tem- 
peratures, the interplay of disorder and interactions in- 
evitably plays an important role in transport. Efforts to 
investigate conductivity near magnetic instabilities often 
rely on the spin-fermion model [3]: the conductivity is de- 
termined by low-energy fermionic excitations, interacting 
with collective bosonic spin modes that carry no charge 
yet become soft modes at the magnetic critical point. 

If the interaction corrections to the conductivity are 
small, they can be analysed in the framework of the spin- 
fermion model using perturbation theory. For instance, 
the interaction contributions to the conductivity of a 2D 
metal near a ferromagnetic (FM) instability can be found 
microscopically [4 using the diagrammatic technique sim- 
ilar to the electron-electron interaction corrections [5j |6] 
to the conductivity of a disordered metal. A separate 
analysis is needed in the case of an antiferromagnetic 
(AFM) instability, characteristic of pnictide supercon- 
ductors, certain heavy-fermion materials and possibly the 
cuprate systems and organic charge transfer salts. Near 
the AFM transition, the momenta of the lowest-energy 
spin fluctuations are large and close to the reciprocal vec- 
tors of the spin superlattice in the AFM phase. As a 
result, electrons strongly interact with spin fluctuations 
only close to narrow regions of the Fermi surface, the so- 
called "hot spots", and can be scattered from one such 
region to another. 

In a system without disorder quasiparticle lifetime and 
weight vanish at the hot spots at low energies. The anal- 
ysis of the Boltzmann transport equation concluded [7], 
however, that the conductivity is dominated by weakly 
scattering "cold spots" . The role of impurities within the 
kinetic-equation approach was investigated by Ueda[8 
and Rosch[9] who found a correction Ap ex T 3 / 2 to the 



residual resistivity po caused by impurities. The frac- 
tional power here was explained by the width of the "hot 
lines" cx T 1 / 2 [9]. It is however unclear why such a 
quasi-classical Boltzmann approach should be applicable, 
considering, in particular, the singular scattering rate 
and vanishing quasi-particle weight. Another interesting 
question, that deserves a separate investigation, is related 
to the notion of impurity-induced " local criticality" [10| : 
the single-particle self-energy may become momentum 
independent E (k,zcj) ~ E (iuj) and show singular be- 
haviour on the entire Fermi surface away from the hot 
spot, that would manifest itself in the resistivity p(T). 
As estimated in Ref. [TT] , such local scattering processes 
lead to the quasiparticle scattering time r s ~ x (T) cx T 3 / 2 . 
This momentum-independent scattering rate implies a 
similar temperature dependence of Ap. The emergence 
of an impurities-induced local single-particle scattering 
rate is particularly interesting given the experimental in- 
dications in favour of local, i.e. momentum independent, 
criticality of electrons near the antiferromagnetic quan- 
tum critical points [2]. 

In this Letter we study microscopically transport in the 
spin-fermion model in presence of impurities using the 
diagrammatic technique. We demonstrate that at suffi- 
ciently strong disorder and weak spin-fermion coupling 
the interaction corrections to the resistivity are domi- 
nated by the processes near the hot spots. The quasipar- 
ticle self-energies in the "cold regions" also show critical 
behaviour, that leads, however, to a smaller contribution 
to the resistivity correction. We find the conductivity 
dependency on temperature and frequency. In the zero- 
frequency low-temperature limit we recover the results 
previously known from the kinetic-equation approach. 

Model. We consider for simplicity a spherical Fermi 
surface. A fermion can absorb or emit a boson only if 
the momentum of the latter is close to a certain large 
vector Q (|Q| ~ /c^), determined by the geometry of the 
AFM spin superlattice in the AFM phase. In order to 
address the experimentally relevant case of an isotropic 
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conduction we assume that there are three such (incom- 
mensurate) vectors, Qi, Q2, and Q3, of equal length 
and directed respectively along different coordinate axes. 
The fermions strongly interact with bosons only close to 
the "hot-spots", points on the Fermi surface separated 
by the vectors Q n , Fig. [I] In our case the hot spots are 
three pairs of circles; an electron near each circle can be 
inelastically scattered to the corresponding circle in the 
opposite hemisphere. 
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FIG. 1: (Colour online) The Fermi surface. 
The action of the spin-fermion model reads 
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where c\ (x) and c a (x) are the Grassman fields for an 
electron (fermion) with spin cr, x = (r,r); S n (x) is the 
bosonic field of the collective spin excitations; n = 1, 2, 3 
labels a pair of hot spot circles separated by the vector 
Q n , g is the coupling constant between the bosons and 
the fermions, £k is the electron spectrum, D n (x — x') 
is the propagator of the bosonic modes; the disorder 
is represented by the random potential J7(r), that acts 
on the fermions only. We use the abbreviation J • • • = 

J d 3 r dr • • • and set e = H = 1. Also, below we sup- 
press the bosonic index n, if the respective expression 
does not depend on it, and use conventions k = (k, iui) 

The energies of all excitations in the spin-fermion 
model are limited by some energy scale A, which is sig- 
nificantly exceeded by the microscopic bandwidth W of 
the tight-binding Hamiltonian of the underlying lattice. 
The shape of the Fermi surface and the excitation spec- 
tra are assumed renormalized upon having integrated out 
all the modes with higher energies A < e < W [3]. The 
relative smallness of the cutoff A allows one to linearise 
the fermionic spectrum £k with respect to k — kp. 

We consider a system sufficiently close to the criti- 
cal point, so that the collective spin excitations are the 
most important bosonic excitations in the model, so one 
can disregard the other types of interaction. On the 



other hand, to address small interaction corrections to 
the conductivity of a disordered metal, observed in exper- 
iments, we assume that the dimensionless coupling con- 
stant, a c± g 2 /vp <C 1, between the fermions and bosons 
is the very smallest parameter in the theory, which allows 
us to treat the interactions perturbatively. In principle, 
the renormalized bosonic propagator D n (iQ,q) depends 
on the strength of disorder and the coupling g. How- 
ever, we assume in what follows that the boson dynamics 
is characterised by phenomenological energy scales inde- 
pendent of g and the elastic scattering time r. 

The assumption of small a implies, in particular, that 
the conductivity is dominated by disorder scattering. 
For simplicity, the impurity scattering is assumed to be 
isotropic, and the disorder potential - weak and Gaus- 
sian; (U(r)U(r>)) = (27rp F r)-^(r-rO, where p F = ^ 
is the density of states on the Fermi surface. The param- 
eter (sft) -1 <C 1 is assumed to be the second smallest 
in the problem, ep being the Fermi energy. 

Perturbation theory. Under the assumptions that 
the bosonic modes are correlated on a short scale yet 
the spin-fermion coupling is weak, one can conveniently 
calculate the interaction corrections to the conductiv- 
ity perturbatively. To the th order in the interac- 
tion, the conductivity is given by the Drude contribu- 
tion <Jo (lj) = (2/3)v 2 pp (t -1 — iuS) , together with the 
weak-localization and the electron-electron interaction 
corrections to it. Next we calculate the leading correction 
(J2 oc g 2 in the spin-fermion coupling g. 

Let us consider first the renormalization of the fermion- 
boson interaction vertex by disorder. The first non- 
vanishing correction to the coupling g, Fig. [2j reads 



5g(q, iu,ie) 
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dk 
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where iuj and q ~ Q are respectively the frequency 
and the momentum of the spin excitation at the ver- 
tex, G(ie, k) = [is — £k + V(^ r ) signs] -1 is the disorder- 
averaged fermion propagator [12]. Because there are no 
excitations with energies larger than A, the integration 
is confined to narrow regions around the hot-spot circles, 
such that |£k,£k+q| < A. This allows one to integrate 
separately the two propagators near the two different hot- 
spot circles, arriving at a very small renormalization of 
the coupling Sg/g ~ (£FT) -1 sign (cj) sign (uo + e), which, 
as we show below, may be neglected when calculating 
the conductivity. Due to the large momentum transfer 
by the AFM bosonic fluctuations, |q| ~ |Q| ~ hp, the 
renormalization is fully perturbative, i.e. may be consid- 
ered only to the leading order in the disorder potential, 
and does not involve the summation of the whole diffu- 
sion ladder, unlike the case of electron-electron or FM 
interactions, when the momentum scattering at the ver- 
tex is small [5] E]. 

Let us notice also that the Hartree contributions to 
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the conductivity and to the quasiparticle scattering rates, 
that contain one boson propagator, (e.g., in Fig.J^J}) van- 
ish due to the symmetry of the spin-fermion interaction 
vertex. Indeed, such diagrams contain an independent 
summation of one interaction vertex with respect to the 
spin polarizations oc ^2 a a a(T ' = 0. 

a) b) 

FIG. 2: a)Renormalization of the spin-fermion ineraction ver- 
tex. b)Hartree contributions to the conductivity and to the 
fermion scattering rate. 
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FIG. 3: Interaction corrections to the conductivity. Diagrams 
1-5 mimic processes close to the hot spots, while 6 is the 
contribution of electrons on the whole Fermi surface. 

To the lowest order in interactions and in disorder 
strength the conductivity correction is given by the di- 
agrams 1-3 in Fig. [3] (thin solid lines correspond to the 
disorder- averaged electron propagators). Their contribu- 
tion to the conductivity is given by the analytic continua- 
tion from the Matsubara to real frequencies, iuo —> u; + i0, 
of the quantity 
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where q = (q, iO) and v\\ if the component of velocity 
at the hot spot parallel to the respective vector Q. The 
analytic continuation reads 



Sa(uj) = 
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where D R (e,q) is the retarded boson propagator, 
A {uj,s) = (uj — e) [n(e) — n(e — uj)] with the Bose dis- 
tribution function n(e). 

Diagrams 4-6 in Fig. [3] contain one extra impurity line 
and represent the next-order corrections to the conduc- 
tivity in the disorder strength, which are not accounted 
for by the renormalization of the interaction vertex. We 
show in what immediately follows that, because of the 
large bosonic momenta, adding more impurity lines to 
the diagrams 1-3 results in small corrections, that are 
perturbative in the disorder strength. 

The straightforward evaluation of the diagram 4 
shows that the corresponding contribution Sa^iuj) ~ 



Sa(iuj)G [(sft) -1 ] is suppressed by the small parame- 
ter (sft) -1 compared to the diagrams 1-3, which can be 
understood as follows. The extra impurity line in dia- 
gram 4 adds to the respective integral one more momen- 
tum integration and two propagators, whose momenta 
are shifted by a large constant vector of the order of kp 
with respect to each other. This extra integration results 
in the relative smallness (e^r) -1 of the diagram. A sim- 
ilar argument proves the smallness of the diagrams with 
crossing impurity lines in the usual disorder averaging- 
diagrammatic technique [12]. Let us emphasize that this 
smallness is specific of the case of the large momen- 
tum scattering by the bosonic modes close to the AFM 
criticality. If the electron dynamics is strongly affected 
by collective excitations with small momenta, e.g., FM 
fluctuations^ or electron-electron interactions [5], then 
taking into account the whole diffusion ladder is neces- 
sary in place of the impurity lines for the diagrams 4 and 
5, as well as in the renormalization of the spin-fermion 
interaction vertex. 

Similarly, the contribution of the diagram 5, 



5a 5 (iuj) ~ ar" 1 ^" 1 + uj) 2 / D(iQ,q) 



(5) 



is also suppressed by the small parameter (e^r) -1 , com- 
pared to the diagrams 1-3. However, Eq. ([5| involves 
the summation of the boson propagator over all Mat- 
subara frequencies, unlike Eq. (J3|, where the frequency 
summation is restricted to \Q\ < \uj\. This difference can 
make diagram 5 important at sufficiently low frequencies 
uj. However, the temperature dependencies of both dia- 
grams are determined by a few terms with Q ~ T in the 
sums in Eqs. ([5| and (J3|, and we can neglect the tempera- 
ture dependency contribution of the diagram 5 so long as 
Eft ^> 1. Sa^(uj = 0,T = 0) amounts to a contribution 
to the residual conductivity due to the spin fluctuations 
and may be disregarded in what follows. Similarly, one 
can neglect the diagrams for the conductivity that cor- 
respond to the renormalisation of the interaction vertex 
since SpT ^> 1. 

Let us proceed to the diagram 6 in Fig. [3j This in- 
teraction correction comes from electrons near the whole 
Fermi surface, unlike the previously considered diagrams 
1-5, that correspond to processes near the hot spots. In- 
deed, anywhere on the Fermi surface an electron can be 
scattered into a hot spot, where its dynamics is impeded 
by the spin fluctuations, cf. diagram 3 in Fig. [4| The 
value of the diagram 6 in Fig. [3] 



5(jQ(iuj) ~ (Ar) 1 5a(iuj) 



(6) 



is small compared to the value Sa(iu) of the diagrams 
1-3 due to the small parameter (Ar) -1 <C 1, which can 
be understood as follows. On one hand, extra smallness, 
~ (£_ft) _1 , comes from the extra impurity line. On the 
other hand, the diagram contains an extra largeness in 
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it, ~ £p/A,- the ratio of the area of the Fermi surface to 
the characteristic size of the hot spots. The combination 
of the two factors leads to Eq. (J6|. 



FIG. 4: Interaction corrections to the scattering rate. Dia- 
grams 1 and 2 describe scattering close to the hot spots, while 
3 is an isotropic contribution to the scattering on the whole 
Fermi surface. 

Thus, we have shown that to the leading order in disor- 
der and interactions the conductivity is given by Eq. (|3| , 
corresponding to the diagrams 1-3 in Fig. [3j To the first 
order in a the conductivity reads 



a(uj) 



2 2 r 

3^ PF [7 



■ io;(l + A)]" 



where we introduced 
6na v 



St- 1 = 



A 



^ J ImD R (e,q)A(w,s) 



ded<\ 

k F uj v J v-'"*'"v~'-'(2^)4' 
6 ™ V -± J ReD R (s,q)A(LO,s) d£dq 



(7) 

(8) 
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k F w 2 v J v v ' y (27r) 4 ' 

The limit of low frequencies and temperatures. To 
make further progress at sufficiently small to and T we 
consider the most general form of the overdamped bo- 
son propagator that follows from the effective Ginzburg- 
Landau-type description: 

D R (s, q) = - [r 2 + (q ± Q) 2 " fefcf.7" 1 ] _1 • (10) 

± 

at \e\ <C e^7 _1 . For \e\ > £^7 _1 , £: " 2 - Here 7 is 

the phenomenological scale that characterises elastic and 
inelastic scattering of bosons. 

Substituting Eq. ([To| into Q yields Sr = 



3aC(3/2) v\\ T 3 / 2 
4V^ 



for (fc F f ) > 7T- 1 and = 0. While at 

(/VO 2 < 7 T_1 and uj = it holds St = 7 f^ T ^f A - 
Similarly, for T = and finite frequencies 5r = 
— ^T72- The crossover between these regimes is de- 
scribed by the expression 
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with the scaling function 

x(z — 
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$(2/,z) = 
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dxdt. 
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In particular, in the dc limit the correction to the resis- 
tivity due to the spin fluctuations at the AFM critical 
point reads 



Sp — aC- 
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)2 k 2 F r 2 ^ 
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where C = 97r — 23.14.., which up to a numerical 
coefficient matches the result obtained within the kinetic- 
equation approach. 

High frequencies. At larger uj the frequency depen- 
dency of the conductivity can be found regardless of the 
particular form of the boson propagator. If the frequency 
is very large, the cutoff A of the excitation energies in 
the spin-fermion model should be chosen larger than uj 
to ensure the existence of quasiparticles that can absorb 
a quantum uj. 

Let us assume first that the bosonic dynamics is deter- 
mined by certain energy scales independent of the cutoff 
A. Then one can neglect the bosonic frequencies Q in 
comparison with the large uj in Eq. (|3j). Making the an- 
alytic continuation and introducing the average value of 
the spin fluctuations (S 2 ) = — J q), we arrive at 



5a(uj) 



avukf 
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The result (14) has a simple physical interpretation. 



At very high frequencies the spin fluctuations may be 
considered frozen and equivalent to static disorder. The 
appropriate modification of the elastic scattering time, 
averaged with respect to angles, estimates S(l/r) ~ 
a(S 2 )/kF and causes correction 5o~ ~ (5r~ 1 )dcro/dr~ 1 
to the Drude conductivity o~q. As necessary, this correc- 



tion matches Eq. (14). 



Intermediate frequencies. Let us proceed to the inter- 
mediate frequencies; uj exceeds the temperature T, but 
is smaller than the smallest characteristic energy of the 
boson dynamics, e.g., the bosonic scattering rate e 2 F ^~ x . 

Then at |fi| < uj the integral J D(ifl, q) ^ 2 °^ 3 is nearly 
independent of because the value of ft affects the bo- 
son propagator only on a sufficiently small momentum 
interval, while the rest of the integral is accumulated on 
a greater interval. 

Thus, in Eq. ^ we may set J D(iQ,q)(dq) ~ 
J D(0, q)(dq). Then the frequency summation in Eq. ^ 
and the analytic continuation to real frequencies yield 



Sa(uj) oc —iuj(r 1 — iuj)~ 



(15) 



independently of the form of the bosonic propagator. 

Discussion. Spin-fermion interactions modify the 
quasiparticle self-energy part on the whole Fermi surface. 
Away from the hot spots the modification 



\Q\)D(iQ, q) signs. 



(16) 

is given by diagram 3 in Fig. [4j It is momentum- 
independent and "locally critical", corresponding to the 
scattering rate r^} oc const + T 3 / 2 . Let us notice that 
Eqs. ([8]), ([9]), and (16) hold for an arbitrary boson prop- 
agator and can be used to analyse transport in, e.g., a 
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system with a two-dimensional spin dynamics [2]. In ar- 
bitrary dimensions d holds r^ 1 oc const + T d l 2 . 

We demonastrated, however, that the interaction cor- 
rection to the conductivity near the AFM instability is 
dominated by the processes near the hot spots, corre- 
sponding to diagrams 1-3 in Fig. [3] We have found the 
dependency of conductivity on frequency uo and temper- 
ature T. In the limit uo = we recover the temperature 
dependencies of the interaction correction to conductiv- 
ity previously known from the kinetic equation analysis; 
5o~ oc T 3 / 2 and 5o~ oc T 2 at the critical point and away 
from it respectively. At T = we find da oc oj 3 / 2 . At 
sufficiently high frequencies the correction is independent 
of temperature and a particular form of the spin propa- 
gator. At very high and moderate frequencies the depen- 



dency is given by Eqs. (14) and (15) respectively. 
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